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Abstract 


A Vlasov  equilibrium  model  is  developed  which  can  describe 
a broad  class  of  beam-plasma  systems  in  an  applied  field,  B^. 

The  model  treats  arbitrary  density  profiles  and  azimuthal 
self-magnetic  fields  as  well  as  a thermal  spread  in  the  momentum 
of  the  beam.  It  is  assumed  that  the  electron  gyroradius  is 
small  but  not  negligible  compared  with  the  beam  radius,  and 
that  the  beam  is  completely  charge  neutralized  but  only 
fractionally  current  neutralized.  The  electron  orbit  equations 
are  solved  and  are  used  to  simplify  the  form  of  the  distribution 
function  and  the  calculation  of  its  velocity  moments.  The 
rotation  frequency  of  the  beam  and  the  diamagnetic  or  para- 
magnetic properties  of  such  warm  beam-plasma  systems  are  also 
studied. 
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I.  INTRODUCTION 


In  recent  years  there  has  been  much  interest  in  the. equilibrium 

and  stability  properties  of  relativistic  electron  beams.  One  particular 

reason  for  this  interest  is  the  use  of  relativistic  electron  beams  as 

l— 8 

a plasma  heating  device.  An  understanding  of  the  physical  mechanisms 

involved  in  this  heating  process  requires  a thorough  knowledge  of  the 

equilibrium  and  stability  properties  of  the  beam-plasma  system.  In 

general  the  stability  analyses  of  such  beam-plasma  systems  that  have 

appeared  in  the  literature  have  assumed  that  the  background  plasma 

6~13 

provides  both  complete  charge  and  current  neutralization. 

For  the  stability  analysis  it  is  also  generally  assumed  that  the  beam  is 
cold  as  well  as  spatially  uniform  and  infinite  in  radial  extent. 

Some  work  has  been  done  where  only  partial  charge  and  current  neutraliza- 

14  1 

tion  have  been  assumed  and  where  finite-radius  beams  have  been  studied 

but  a more  complete  analysis  is  needed.  The  presence  of  self  magnetic 

fields  and  radial  boundaries  can  have  Important  effects  on  the  mode 

structure  of  beam-plasma  systems;  for  example,  their  presence  can 

greatly  enhance  the  coupling  between  the  transverse  and  longitudinal 

waves  that  are  found  by  treating  a uniform,  current  neutralized  beam. 

The  purpose  of  this  article  is  to  present  a Vlasov  equilibrium  model 

for  a warm  relativistic  electron  beam  propagating  through  a background 

plasma  immersed  in  a uniform  applied  magnetic  field,  B^.  The  model 

allows  for  self-magnetic  ffelds,'  for  arbitrary  radial  beam  density 

profiles  and  for  a thermal  spread  in  momentum,  however,  it  will  be 

assumed  that  the  beam  is  charge  neutralized  by  the  background  ions. 

Since  Bg^O  and  the  beam  is  warm,  both  paramagnetic  and  diamagnetic  beam 

equilibria  are  possible,  where  the  magnitude  of  B0  can  be  adjusted  by 


varying  the  fractional  current  neutralization  provided  by  the  counter- 
streaming  of  the  background  electrons.  In  a subsequent  paper  this 
equilibrium  model  will  be  used  in  analyzing  the  stability  of  such  a 
beam-plasma  system.  It  should  be  pointed  out  that  this  model  can  easily 
be  generalized  to  describe  the  problem  of  a Tokamak  plasma  with  a supra- 
thermal  electron  population.  Here  again  both  the  toroidal  and  the 
poloidal  magnetic  fields  can  be  described  as  well  as  the  finite  radial 
density  profile. 

The  applied  field,  Bz>  is  assumed  large  enough  so  that  the  gyro- 

radius  of  both  the  beam  electrons  and  the  background  electrons  is  small 

compared  to  the  radius  of  the  beam.  With  this  assumption  the  electron 

orbits  are  easily  calculated  by  an  expansion  method  using  the  smallness 

20 

of  the  gyroradius  and  the  method  of  harmonic  balance.  These  orbits 

are  useful  in  writing  down  the  distribution  functions  for  the  two 

electron  components  in  a convenient  form  that  will  facilitate  the 

calculation  of  the  velocity  moments.  The  ions  will  be  described  by 

a stationary  cold  fluid  model  which  will  be  shown  to  be  appropriate 

m or  r? 

for  a>>>u)  . and  T.  < £ . 

ci  i ^ 2 

In  Section  II  a discussion  of  the  assumptions  that  are  appropriate 
in  describing  this  beam-plasma  system  will  be  presented.  Following  this 
discussion,  the  procedure  for  calculating  the  electron  orbits  will  be 
outlined  in  Section  III.  The  beam  electron  distribution  function,  f,  , 

D 

and  its  velocity  moments  are  developed  in  Sections  IV  and  V,  whereas 
Sections  VI  and  VII  provide  a similar  development  for  the  background 
electrons.  The  equations  that  describe  the  self-consistent  magnetic 
fields  are  presented  in  Section  VIII  and  the  numerical  solution  of  these 
equations  are  included  in  Section  IX,  which  gives  some  examples  of 
typical  beam-plasma  equilibria  described  by  this  model.  Finally, 
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vn  B0 

lines,  — , and  the  diamagnetic  drift.  By  requiring  the  ExB  drift  to 

B 

z 

be  small  in  comparison  to  the  azimuthal  streaming  the  following  condition 


on  the  ion  temperature  is  obtained 


T.  <<  eB.r,  — 
i 6 b V1^ 


Here  v(|  /c,'0(l)  was  assumed  for  the  beam  electrons.  If,  for  example, 

BgM.  kg,  r^l  cm  and  nj/nj)'v'25.  Eq.  (6)  yields  ^<<7.5  keV,  which 

is  easily  satisfied  in  typical  beam-plasma  systems.  In  addition, 

if  the  time  scale  of  interest  in  a stability  analysis  is  much  shorter 

than  u>~ (i.e.,  w>>0)  J) , then  the  ions  can  be  considered  as  cold  if 
ci  ci 


v , < a)  .r.  <<  u»r.  . 
th  'v  ci  b b 


This  can  also  be  written  as  a condition  on  as 

2 2 2 2 
m.o)  r.  , m a)  r. 

^-M^"|keV  (8) 

where  Eq.  (2)  has  been  used.  Note  that  this  is  in  agreement  with  the 
example  given  for  Eq.  (6).  Thus  for  typical  conditions  of  beam-plasma 
experiments  it  is  reasonable  to  treat  the  ions  as  a cold  fluid  for 
the  purposes  of  studying  the  stability  of  waves  with  o)>>wc^.  This 
frequency  condition  also  implies  that  the  ions  are  unmagnetized. 

Since  J£=0  the  relativistic  factor,  y,  for  the  electrons  is  a constant 
of  the  motion.  This  fact  along  with  the  assumption  that 


where  p&  is  the  electron  gyroradius,  simplify  the  calculation  of  the 
electron  orbits.  These  orbits  will  be  invalid  for  a small  fraction 
of  the  electrons  within  a few  gyroradii  of  the  axis  but  will  properly 


describe  the  bulk  of  the  electrons.  It  should  be  noted  that  the  beam 


density  gradient  and  the  self-magnetic  field  in  the  0 direction  both 
are  large  at  the  edge  of  the  beam  where  these  orbits  are  appropriate. 
Therefore,  the  effects  of  the  self  field,  B_,  and  the  finite  radial 


size  of  the  beam  can  be  studied  in  detail. 


III.  ELECTRON  ORBITS 


The  calculation  of  the  electron  orbits  is  greatly  simplified 
when  the  electric  field  is  negligible  and  the  electron  gyroradius  is 
much  smaller  than  the  beam  radius.  Under  these  conditions  the  equation 
of  motion  can  be  written  as 

y di = "(r^ce)  (10) 

2 eS 

where  H“Ym  c is  a constant  of  the  motion  and  m = . In  addition 

e ''ce  m c 

e 

to  H,  the  canonical  angular  momentum,  PQ , and  the  canonical  z-momentum, 

P , are  also  constants  of  the  motion  since  it  is  assumed  that  there  is 
z 

no  functional  dependence  on  the  6 or  z coordinates.  Solving  Pg  for  0 


and  P for  z one  obtains 
z 


de  s pe  . “e 

T - 1,9 2+~ 

m r 
e 


— + a 
m z 
e 


where  a.  = — and  A.  is  the  j component  of  the  vector  potential 
3 m c j 

• • 

(j«0,z) . These  expressions  for  0 and  z are  then  used  in  the  r component 


of  Eq.  (10)  to  obtain  an  equation  for  r(t): 


2- 

Y r - r 


+ a J a)6  -(  — + a ) o)Z  = F(r)  (13) 

z / ce  \ mr  0 / ce 


If  r(t)  is  written  as 


r(t)  - R(l+p(t))  , 
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then  F(r)  in  Eq.  (13)  can  be  Taylor  expanded  about  r=R  using  the  assump- 
tion that  | p | <1 , where  | p | R is  Identified  with  the  electron  gyroradius 
and  R is  identified  with  the  electron  guiding  center  position.  Again 
it  should  be  stressed  that  these  orbits  do  not  correctly  describe 
the  motion  of  a small  fraction  of  the  electrons  near  the  axis  of  the 
beam;  there  | p | >1 , so  that  the  Taylor  expansion  would  not  be  appropriate. 
This  is  not  a serious  problem  as  long  as  r^>> | p | R,  so  that  the  orbits 
do  properly  describe  the  bulk  of  the  electrons.  Expanding  F(r) , Eq. 

(13)  becomes 


Y p = <Xq  + a^p  + o^p  + . . . 


whete 


00  n n- 
y fL_  i_Z 

n=0  n'  ar11  r=R 


2 

X„  X„  X 

0 . 0 Z Z 0 

Cln  ■ o + 0) 0) 

0 L r ce  r ce 


0 „ 6 


c Xj  * 5a 


o 0 , 3u  / 3m  3 a)  \| 

a,  * (ai  ) + (a)  ) - 3a~  + x„  — x l + 1 

1 ce  ce  0 0 3r  z \ r 3r  / r=R 

2  z z (17) 

. *e  U1  8 “ce  1L  8“ce  . e \ 

>0  + T (T  - 2r  W * “ce) 

(»‘  )Z+Tr  «•'  >2  + (»*  )2)  '~A(!r  ■ 5~e  )1 

ce  4 3r  ce  ce  r \ 2 jr  ce'l 


ce  _ ce  . z 
~2  2r  3r  + “ce 


• « p • 

x0  .^  + VR> 

e 


x = — + a (R)  . 
z tn  z 
e 


Note  that  the  relationship 
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has  been  used  in  deriving  the  expressions  for  a^,  and  a^. 

Since  F(p=0)=agyo,  Eq.  (15)  allows  for  corrections  to  the  guiding 
center  position,  R.  Unlike  traditional  guiding  center  theory  which  does 
not  treat  the  problem  of  streaming  electrons,  this  correction  to  R 
can  be  as  large  as  0(e)  for  the  case  where  v((  /c'UD(l).  The  correction 
can  not  become  of  0(1)  since  it  was  assumed  |p|<l.  If  a generalized 
potential,  V,  is  defined  through 

F(p)  " - U • <22> 


it  is  easy  to  see  that  this  correction  to  the  guiding  center  position 

can  be  attributed  in  part  to  a small  asymmetry  in  the  potential  and 

in  part  to  a small  shift  away  from  p=0  in  the  minimum  of  the  potential 
22 

well. 

For  convenience  the  Taylor  expansion  in  Eq.  (15)  will  be  terminated 

after  the  quadratic  term.  The  resulting  nonlinear  equation  can  be 

23 

solved  analytically  in  terms  of  elliptic  integrals,  however  that 
procedure  does  not  provide  any  physical  insight  into  the  solution. 
Insight  can  be  obtained  by  writing  p as 


where 


p = pQ+p1sin4>+P2Cos2(J>+. . . 
<j>=fit+(f>Q  . 


(23) 

(24) 


Here  PqM)(e),  p^0(e)  and  P2'''0(e  ),  where  e<<l  and  is  used  to  order  the 
variables.  Substituting  Eq.  (23)  into  Eq.  (15)  one  obtains 


/ 2 n2pl\  2 2 

V*0  + ot^  + o^Pq  + — j—)  + (-8  Y +a1+2a2P0)p]Lsin<J> 

2 

/ 2 2 a2p  1 \ 1 
+ (j-4n  y p2  + 0^2 2 — / C°s2<j>-K)(e  )=0 


(25) 


2 


J 


This  equation  can  then  be  solved  using  the  method  of  harmonic  balance. 
In  order  to  further  simplify  the  solution,  each  variable,  s,  will  be 
written  as 

(j)  (26) 


s = I s 
j 


.<j> 


where  s''J',  will  denote  the  contribution  to  s which  is  of  0(r'*). 

2 

In  this  manner  the  solution  of  Eq.  (25)  through  0(e  ) is  reduced 
to  solving  the  following  set  of  equations: 


(°)  n 
a0  0 


♦ w - » 


(0) 


(p,(1))2  = 


-YV)(0)+«<0)=0 


-YV)(1)  +a<»  +2a(2°V11)'0 

JO) 


-4y2(Q2)(0)p<2)  +4°V22> 


V (P^)2  = 0 


(27) 

(28) 

0 (29) 

(30) 

(31) 

(32) 


where  Pq^*p^=p  by  assumption. 


one 


Ordering  the  variables  appropriately  to  describe  the  beam  electrons, 

PI1  P±  XZ  X0  1 

has  'vO(l),  ^O(c),  — ^0(1)  and  — ^0(e;  along  with  the 

me  me  ’ c c 

e e p, 


to 


assumption  of  Eqs.  (1)  and  (2).  Then  setting  4> ( t=0)  = 4>q  and  P^=  ^ 

e ' 

satisfy  the  initial  conditions,  Eq.  (15)  is  solved  for  p(t).  It 
should  be  mentioned  that  the  correction  to  the  guiding  center  position, 
R,  is  0(e)  for  the  ordering  considered  here  (i.e.,  Pq^^O).  This 
point  is  discussed  in  more  detail  further  on  in  this  section.  By  using 
this  expression  for  p(t),  the  momentum  of  the  beam  electrons  can  be 


obtained 
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Pr  = vx 


cos  9 - 


p9  = 


D OJ 

ce 


_ z 

2m  rco 
e ce 

2 

i P. 


sin2 -f-  0(  e ) 


2 


Pm  («;j 


62 
ce 


ce 


+ p,  ^1  4 
2 


2 z 

m rto 
e ce 


P i w 
K 11  ce 


, z . 3 
m r (to  ) 
e ce 


ll 


3to 


ce 


2m  to2 
e ce 


2m  r(toZ  )2  2m  (toZ  )2  9r 
e ce  e ce 


cos2<f>  + 0(e  ) 


sintf) 


pz  = p il 


with  <t>  = fit+4>Q  = 


1+J 


ce 

Y 


0 2- 
o 

ce 

z 
o 

ce 

0 

p ,,  to 
v 11  ce 

2 

2m  rvio 
e ce 


Piw  o 

— sin<)>  + 0(e  ) 


ce 


3to 


3 + 


ce 


0 3r 
o 

ce 


t + 


(33) 


(34) 


(35) 


(36) 


Here  II  refers  to  the  streaming  (<{>  independent)  component  of  the  motion, 
l refers  to  the  gyromotion,  and  x^^  has  been  identified  as  p |(  /m^y. 

Eq.  (15)  was  solved  with  all  the  variables  written  as  functions  of  R; 
then  the  intermediate  result  for  j)  was  expanded  about  R=r  in  order  to 
arrive  at  the  solutions  exhibited  in  Eqs.  (33),  (34)  and  (35),  where  all 
the  variables  ar|  again  written  as  function  of  r.  Also  note  that 


Q is  not  simply 


ce 


but  contains  a drift  correction  of  0(e), 


From  Eqs.  (34)  and  (35)  it  is  clear  that  the  streaming  component 
of  the  motion  involves  both  streaming  along  the  field  lines  as  well 
as  a drift  across  the  field  lines.  Defining  a pitch  angle,  n,  for 


the  streaming  motion  one  has 


Pl(  sinn 


e 

PH  “ce  ^ 1 
z 2 


2 

P, 


2.6.2 
p,i  (“ce) 


ce 


m rwZ  m r(wZ  )^ 
e ce  e ce 


p(j  cosq  = p( 


/ 0 2' 

( ,/(!)' 

1 +i(— e 

V 2\(.z 

ce 


(37) 


(38) 


so  that 
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0) 

ce 

z 

0) 

ce 

' 2 

Pi 

_i_ 

D2  1 

/ 0 
/ U) 

' ce 

2 + 

Pll  ' 

[ Z 
ce 

m p,.  ru  \to  / I 

e H ce  ce  J 

Thus  the  streaming  is  predominantly  along  the  field  lines  for  the 

case  considered  here,  however  the  streaming  does  depart  slightly  from 

(2) 

this  course  due  to  the  curvature  drift  as  indicated  by  n 

The  fact  that  the  0(e)  correction  to  the  guiding  center  position, 

Pq1^ , is  nonzero  for  this  ordering  can  be  explained  by  consideration 

of  the  conservation  of  canonical  angular  momentum,  Pn.  If  each  beam  elec- 

tron  is  injected  axially  into  the  plasma  at  some  r=R  with  zero  mechanical 

angular  momentum  from  a source  which  is  immersed  in  the  axial  applied 

magnetic  field,  then  the  canonical  angular  momentum  of  the  electron  is 
z 2 

simply  -tl)ceR  /2.  The  contribution  of  the  axial  self-magnetic  field  to  PQ 

is  not  considered  here  since  it  is  two  orders  of  magnitude  smaller  than 

the  contribution  of  the  applied  field.  Now  as  the  electron  begins 

streaming  along  the  helical  field  lines  while  gyrating  about  them,  it 

gains  mechanical  angular  momentum  at  the  expense  of  the  angular  momentum 

associated  with  the  applied  field.  In  terms  of  anharmonic  oscillation 

potential,  Eq.  (22),  the  potential  minimum,  R,  and  the  phase-averaged 

particle  position,  R,  differ  in  first  order,  by  an  amount  which  depends 

on  the  injection  v , since  the  potential  is  a velocity-dependent  one. 

• - • 

The  rotation  associated  wj.£h  the  streaming  along  the  field  lines  accounts 

for  the  largest  part  of  the  mechanical  angular  momentum  and  is  of  0(c). 

This  process  is  equivalent  to  saying  that  the  guiding  center  of  the 

- 0 z 2 

electron  is  injected  at  R=R+v..  w /(w  ) , and  not  at  R.  The  electron 

J II  ce  ce  ’ 


is  injected  into  the  potential  well  at  R,  displaced  from  the  minimum 


12 


at  R;  thus  the  electron  oscillates  about  R with  a frequency,  ft.  Since  V(p) 
is  not  a simple  harmonic  potential  there  are  also  small  corrections  to 
the  usual  electron  gyrofrequency  as  shown  in  Eq.  (36).  In  what  follows  R 
will  be  used  to  denote  the  electron  guiding  center  position. 

It  is  also  useful  to  write  down  the  single  particle  orbits 
in  the  equilibrium  fields.  These  orbits  are  necessary  in  calculating 
the  perturbed  distribution  function  when  using  the  method  of  charac- 
teristics in  a stability  analysis.  If  at  T=t'-t=0,  the  beam 

electron  orbits  are 


Pr  2 

r ' = r H (sin^-sintfO  + 0(e  ) 

z U 


(41) 


m u) 

e ce 


0 + 


P,,  “ 
r il  ce 


,m  *yxui 
' e ce 


P.L  2 

(cOS<(>-COSi|>q)  + 0(e  ) 


z 

m ru 
e ce 


z ' = z + 


m Y 
e 


T + 0(E) 


(42) 

(43) 


where  terms  of  0(e  ) have  been  dropped  and  where  4>=>ftT+<JiQ . 


In  order  to  describe  the  background  electrons  the  following 

x 

II  / X / 

ordering  is  used: 


O O ^ _ O O 

^0(e  ),  — ^0(e  ),  — ”vO(e  ) and  — ^0(e  ) 

c c c c 


(along  with  the  assumptions  found  in  Eqs.  (1)  and  (2)).  In  this 
case  the  velocity  of  the  electrons  is  found  to  be 


v 

r 


v 

z 


4 

= vAcos<J>  + 0(e  ) 

6 

v„  u . 

= — — Cg  + v,sin<t)  + 0(e  ) 
(0 

ce 

= v ||  + O(e^) 


(44) 

(45) 


(46) 
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where  <t>=<D  t+4> n.  This  shows  that  the  background  electrons,  which  are 

ce  u 

responsible  for  the  partial  current  neutralization,  essentially  stream 
along  the  field  lines  and  gyrate  about  them.  Again  it  is  useful  to 
write  down  the  orbits,  which  for  the  background  electrons  are  simply 
given  by 


1 3 

r ' = r + (sincfi-sinit’  )+0(e  ) 


'0' 


ce 


- 3 

6'  = 0 - (cos<t>-cosi}>n)  + 0(e  ) 


rco 


ce 


z'  = z + v t + 0(e  ) 


(47) 


(48) 


(49) 


where  J^,=^  when  x=0 . 

Although  Eqs.  (27)  through  (32)  have  only  been  solved  here  for 

two  particular  orderings,  other  solutions  are  possible  if  the  variables 

u>0 

are  ordered  differently.  For  example  if  the  order  of  — — is  reduced  to 

2 wce 

0(e  ) instead  of  0(e),  it  is  easy  to  find  the  new  solutions  directly 

from  the  solutions  already  discussed  by  merely  changing  the  ordering 

3 

where  appropriate  and  dropping  any  terms  of  0(e  ) or  smaller.  For 
the  particular  example  mentioned  here,  the  diamagnetic  current  will 
now  be  of  the  same  order  of  magnitude  as  the  azimuthal  current  due 
to  the  electron  streaming  along  the  helical  field  lines,  i.e.. 


0 V 2 

. \ 8£nn  _ 2X 

Pul  ~Z~i  ^ — — — ^ 0(e  ) 


(0  1 

ce 


z 3r 


(50) 


m you 

e ce 


C0 

This  implies  that  the  beam  can  be  diamagnetic  whereas  if  'v,  0(e) 

0) 

ce 

the  beam  is  paramagnetic.  The  diamagnetic  current  is  of  course  a 
fluid  propery  and  does  not  appear  in  the  single  particle  orbits. 


i 
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IV.  THE  BEAM  ELECTRON  DISTRIBUTION  FUNCTION 

Consider  the  distribution  function 


£b  - Kb(R)exp  {-  i IH-4p0-V^P2]}  s Np 


-®b/Tb 


(51) 


Here  w can  be  associated  with  the  fluid  rotation  of  the  electron 
R 


r,b 


beam  and  V can  be  associated  with  the  axial  streaming  of  the  beam, 
z 


Since  the  guiding  center  position  of  the  electron  is  a constant  of 
the  motion,  N can  be  written  as  a function  of  R.  Furthermore, 
if  p](  is  defined  as 


P’.l  = Pll  "5I 


(52) 


such  that  P)(  is  the  average  value  of  the  streaming  component  of  the 


momentum  and  pj.  represents  the  thermal  spread  about  this  average  value, 


then  y can  be  expanded.  Assuming  'v  0(1)  and  ^ ^ 0(e)  one 

me  mcmc 

e e 

finds  that 


2^.  2 /-2 


Y + 


P IJ  P II  . Pi^Il  ^ 


with  Y= 


1 + 


2 2 
m c 
e 


1/2 


-22 

ymec 


o"2  2 

27m  c 
e 


+ 0(eJ) 


(53) 


. Using  this  expression  for  y and  the  expressions 
found  in  Eqs.c(33),  (34)  and  (35)  for  the  components  of  g^  can  be 
written  as 


2,2  .-2 
2-  P-L-^M  ' Y 
gb  = me  y + 


27m 


b 

- mRr 


P„  u 


ce 


+ pisimj)-mea9 


ce 


(54) 


- Vb 
z 


fp .% 

(£*'] 

e 

p.U) 

1 ce 

PU  2 

\ z ) 

z 

\U)  / 

CO 

' ce 

ce 

sin4>  - m a 
e z 


+ 0(e3) 


where 


• 1 


• 9 

...u 


15 


b 

“R 


R;4+5‘H(4))) 

ce  ce 


and  — — ^ 0(e).  By  choosing  N,  (R)  as 
z b 


ce 


•,/,  +i.  (R)/t, 

Nb(R)  = ^(R)(2n^neT)"J/^e  0 


(55) 


(56) 


where 


ib(R)  = 


r 2-  b “ce  \ 

mc  y _ WrR  — - mea6 


ce 


/ 52  0 2 v\ 

(-7s)  - vJ 
(0  / 
ce 


(57) 


r=R 


f,  can  be  written  in  a convenient  form.  To  accomplish  this  g,(R)  and 
b b 

n^(R)  are  Taylor  expanded  about  R=r  where 


P±  2 

R = r sin<{>  + 0(e  ) 

m a) 
e ce 


(58) 


Here  again  it  is  crucial  that  the  gyroradius  be  small  in  comparison 
to  the  radius  of  the  beam.  The  final  result  is 


m a) 
e ce 


nb  , 

z 9—  Sln<fl 


(2irYmeTb) 


-3/2 


exp 


p Win1  ‘ 


(59) 


1 2rme\  I 

where  n (r)  is  an  arbitrary  function  of  r and  is  obviously  identified 

3tlb 

as  the  particle  density.  Moreover,  the  term  involving  sin^i 
accounts  for  the  diamagnetic  current  due  to  the  gyromotion  in  the 
presence  of  a density  gradient. 

The  fact  that  n,(r)  is  an  arbitrary  function  of  r is  of  particular 

D 

interest,  since  it  allows  one  to  model  any  experimentally  determined 
density  profile  which  satisfies  the  condition  that  the  gyroradius  is 
small  in  comparison  with  the  density  gradient  scale  length.  In  Cartesian 
coordinates  an  appropriate  constant  of  the  motion  analogous  to  the  guiding 
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center  position,  R (see  Eq.  (58)),  is 

YV, 


£ - x - 

X 0) 


(60) 


ce 


In  this  geometry  the  appropriate  expansion  for  some  function,  h(£x). 


would  be 


YV.. 


h(£x)  % h(x)  + (-  -^-L) 


3h 


3£ 

ce  x 


(61) 


£ -x 
x 


where  again  the  gyroradius,  | '»'Vy./u>ce  I , must  be  small  in  comparison 


with  the  gradient  scale  length  of  h(£x). 


V.  VELOCITY  MOMENTS  OF  ft 


Using  the  expression  (59)  for  the  beam  electron  distribution 
function,  one  can  easily  determine  the  fluid  properties  of  the  beam. 

For  the  sake  of  convenience,  the  p(|  , pA  and  $ momentum  space  coordinates 


introduced  in  Sec.  Ill  will  be  used  Instead  of  the  p , p.  and  p system. 

r u z 


The  Jacobian  for  this  change  of  variables  is  given  by 


<Vpedpz  “ 


Pi 


l + 


r(" 


dto 


- r 


ce 


2mer4°c/ 

dp1dp,|  d<J> 


3r 


)- 


Pi 


z 

m rw 
e ce 


8in$ 


+ 0(e3) 


(62) 


With  this  change  of  variables  and  with  the  components  of  ^ given  in  Eqs. 


(33),  (34)  and  (35),  the  moments  of  f,  are  readily  calculated.  The 

D 


density,  n^,  the  fluid  velocity,  , the  fluid  momentum,  £b>  and  the 


pressure  tensor,  n^,  are  defined  as  follows: 
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“b". 

3 

fb  d p 

(63) 

^b-. 

f J_  f .3 

meY  b p 

(64) 

nb^b  = j 

* fb  d3p 

(65) 

u 

W (A  - 4,)  fb  d3p  • 

e 

(66) 

The  density  Is  trivially  given  by 

Vr)‘Vr) 

and  the  fluid  velocity  and  momentum  are  found  to  be 


Pbr  = VVbr  = 0 + °<£  > 


(67) 


pbe  “ V\e  " 


p J 

H ce 


th 


ce 

e 


~ z 3r 
2m  u 
e ce 


3u) 


p2  0 

+ ^-3  ( 

, 2 ,3\ 
m r(u>  ) 
e ce 


P.  - m yV,  - PM  + 
bz  e bz  II 


3w 


0 


(68) 


ce  r ”ce)  + o(E3) 


2 3r 


5?. 


0 

0 , z- 72  V^ce 

2mer(u)ce) 


3gj 


- r 


ce 


3r 


) + 0(e2) 

(69) 


where  y is  defined  in  Eq.  (53),  and  p s2m  yT,  . The  reason  for 

th  e b 

3 2 

calculating  P.  to  0(e  ) while  only  calculating  P.  to  0(e  ) is  clear 

Do  DZ 

where  one  examines  the  radial  force  balance  equation,  keeping  in  mind 

<D0 

eg 

that  'v  0(e).  This  equation  is  simply 

z 

0)  o 

ce  p2 

b0  . z „ 0 „ . 


0 — + P,  , - PL  + -J-  [v*n.  ] 

mgr  ce  b0  ce  bz  nb  •v  £b  r 


(70) 


Calculating  the  elements  of  the  pressure  tensor,  results  in  the 
following  expression  for 


« 
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■I 

CL 


W T*r~  '•  T~~r~*Trr~ 


,ee 


^•Vr  = 7l?Kr)  --T 


(71) 


P u 3nu 
*th  b 


2v  9r 


Thus  by  substituting  Eqs.  (68),  (69)  and  (71)  into  Eq.  (70),  one  can 


easily  verify  that  the  radial  force  balance  equation  is  satisfied 
3V 


to  0(e). 

Defining  a pitch  angle,  £,  for  the  fluid  motion  one  has 


Pb0  * Pll  sin^ 


Pb2  = Pll  cos^ 


(72) 

(73) 


where 


P2 


_ = v H / 0 ^ce  \ 

II  ” II  7 . z .2  \wCe  T 3r  ' 

2m  r(w  ) 


P..  = 


- r 
ce 


(74) 


ce 


so  that 


.(1) 


ce 

z 


(75) 


ce 


and 


.(2) 


<"ce>2 


th 


8 £nr 


m PH  r(toz  )^  2m  P..  wZ 
e M ce  e II  ce 


(76) 


To  first  order  the  fluid  also  flows  along  the  field  lines,  as  did 
the  electron '8  guiding  center.  To  second  order,  however,  the  motion  of 
a fluid  element  differs  from  the  trajectory  of  the  guiding  centers. 

This  difference  arises  from  the  electron's  gyromotlon  about  the 
guiding  center  position  and  from  the  presence  of  the  radial  density 
gradient. 


Considering  Eqs.  (75)  and  (76),  the  fluid  motion  of  the  beam  has 

(1) 


a simple  interpretation.  As  mentioned  above  £ shows  that  the  fluid 


predominantly  flows  along  the  helical  magnetic  field  lines.  £ 


(2) 


on 


4 1 
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the  other  hand,  describes  the  shear  of  the  fluid  flow  away  from  the 
field  lines  due  to  curvature  drift  and  the  diamagnetic  drift.  Thus 

3 

the  self  magnetic  field,  B , will  be  paramagnetic  as  long  as  the 
diamagnetic  drift  remains  small.  In  order  to  have  a diamagnetic 
beam  the  diamagnetic  drift  must  be  at  least  as  large  as  this  streaming 
along  JJ,  i.e. , 


P — < 
*11  z £ 

0) 

ce 


th 


9 Inn 


- z 3r 
2m  u 
e ce 


(77) 


This  usually  requires  that  the  current  neutralization  must  be  fairly 
complete.  This  will  be  discussed  in  more  detail  in  Sec.  VIII. 


VI.  THE  BACKGROUND  ELECTRON  DISTRIBUTION  FUNCTION 

The  background  electron  distribution  function,  fe>  is  similar  in 

form  to  f.  given  in  Eq.  (51) 
b 

_ g 

fe  - {-  i i V * ' (78) 

0 

Here  wR  can  be  associated  with  the  fluid  rotation  of  the  background 
electrons  and  V^  can  be  associated  with  axial  counterstreaming 
of  the  background  electrons.  This  counterstreaming  is  responsible 
for  the  partial  current  neutralization.  Since  the  guiding  center 
position  of  the  electron  is  a constant  of  the  motion  Ng  and  V^ 
must  be  written  as  a function  of  R so  that  the  return  current  is 
restricted  to  the  beam  channel  even  though  the  background  electron 
density  varies  on  a length  scale  much  longer  than  r^. 

Using  the  expressions  found  in  Eqs.  (44),  (45)  and  (46)  for  the 
components  of  ^r,  g e becomes 


m 


■7  -v 


M >2«PVe+2Vz-V8  1 + 0<t5) 


(79) 


A R 1 

where  V„  iVe(r)+0(e  ) and  ^ 0(e  ).  By  choosing  N (R)  as 

II  z z c 

-3/2  +ie(R)/Te 
Ne(R)  = ne(2rrmeTe)  e 


where 


ie(R)-2RVe«azV«-v2 


f can  be  written  in  the  following  convenient  form, 
e 

n v^Kv..  -vii 

f • exp  - ; C82) 

e (2rrm  T ) 3/2  v 

2 e e th^ 

with  v L = — - . In  arriving  at  this  result  g (R)  was  Taylor  expanded 
th  m e 

e 

about  R=r  where 


R = r - sin<f>  + 0(e  ) (83) 

z 

w 

ce 

Note  that  ii  may  be  a weak  function  of  r but  is  essentially  constant 
e 

across  the  radius  of  the  beam. 


VII.  VELOCITY  MOMENTS  OF  f 


Using  the  expression  for  the  background  electron  distribution 
function,  ffi,  given  by  Eq.  (82),  the  fluid  properties  of  the  background 
electrons  can  be  investigated.  Again  the  v^  , vA  and  <p  velocity  space 
coordinates  developed  in  Sec.  Ill  will  be  used  in  place  of  v^,  v^  and 
v . The  Jacobian  for  this  change  of  variables  is  simply  given  by 


4 

dvrdv0dvz  * [vi+0(e  ) ]dvAdv(|  d<J> 


With  the  components  of  % in  Eqs.  (44),  (45)  and  (46)  expressed  in 
terms  of  these  natural  coordinates,  the  moments  of  fg  are  readily 
calculated.  Since  the  background  plasma  is  nonrelativistic,  the  density. 
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n , the  fluid  velocity,  ^ » and  the  pressure  tensor,  n , are  defined 


as  follows: 


n = 
e 


n V 
e've 


f d3v 
e 


v f d v 

-v  e 


II  = m (v-V  ) (v-V  ) f d v 

e <1  "ve  'v  ^e  e 


Performing  the  appropriate  integrations  one  finds 


n = n 
e e 


V = 0 + 0(e) 
er 

0 

V..  w / 

Ve0  = + 0(e4) 

0) 

ce 


V - V„  + 0(eH) 
ez  >1 


n = 0 + 0(e) 


(85) 

(86) 
(87) 


(88) 

(89) 

(90) 

(91) 

(92) 


Thus  fluid  motion  of  the  background  electron  only  involves  the  slow 
counterstreaming  along  J3;  all  other  fluid  motions  are  ignorable. 

The  fact  that  V is  an  arbitrary  function  of  r allows  the  freedom  to 
choose  the  return  current  profile  and  thus  determine  the  dependence 
of  B0  on  r. 


VIII.  MAGNETIC  FIELD  EQUATIONS 


The  result  of  Sec.  V and  VII  can  now  be  utilized  to  self  consistently 
determine  the  magnetic  field,  JJ.  From  Maxwell’s  equations,  £ is  given  by 


(93) 


where  is  current  produced  by  the  fluid  motion  of  the  electrons. 
Since  the  ions  are  stationary,  Eq.  (93)  simply  becomes 


'v.  'tee 


2 2 v 

U)  , O) 

2 ^b  2 've 

c c 


where  ai 


4nn  e 

a 


is  the  plasma  frequency  of  species  a.  Substituting 


Eqs.  (68)  and  (69)  for  V,  and  Eqs.  (90)  and  (91)  for  V , the  equations 

'mj  'ue 

for  the  two  components  of  ^ are 

2 _ q 2 

3wZ  as  , (P,.  w • p 3£nn, 

ce  _ pb  ll  ce th  o 

3r  2 z - 2-  z 3r 

c m (o  y 2m  yw 

e ce  e ce 


fn  ("L>2 

2m^Yr(uZ  )^ 
e ce 


r 

3w9  V 

ce  | 

0 

U) 

ce 

)\ 

2 6 
w rV„  u 

i Pe  II  ce 

2 z 

c w 

ce 


13  , e,  pb  rn  . II  ce  /,  r ce  l , pe  r„  , 

717  2 + 7T- ' V ,2  I1  - — 17“  + 2 tvn] 

c (m  v 2tn  Yr(u)  ) \ u /J  c 

e e ce  ce  (96) 

Using  the  assumptions  found  in  Eqs.  (1)  and  (2)  along  with  the 
P 

assumption  that  — - 0(1),  Eq.  (96)  shows  that 

6 2 2 2 

r,  id  , r, 

-^  = -|<0(e)  (97) 

c r 

s 

where  rg  is  skin  depth  of  the  beam.  This  implies  that  the  fields  can 

readily  penetrate  into  the  interior  of  the  beam.  Furthermore,  if 

the  current  neutralization  is  considered  poor,  the  first  term  on  the 

right  hand  side  of  both  Eq.  (95)  and  Eq.  (96)  provides  the  dominant 

WPb  nb 

contribution  to  the  current.  In  this  case  — r—  = — n.  0(e);  however,  if 

z n 

oipg  e 

the  background  plasma  density  is  increased  (or  equivalently  if  the  beam 

V.l  2 

density  is  decreased)  while  — — remains  of  0 (e  ) for  the  background 


electrons,  then  the  current  neutralization  is  improved  although 
it  is  still  not  necessarily  complete.  When  this  occurs  the  situation 
is  more  complicated.  If  the  current  neutralization  is  improved  (e.g.. 


I 
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<50%  complete)  by  increasing  n£,  but  at  q/u  g is  still  0(e),  then  the 


first  and  last  terms  on  the  right  hand  side  of  Eq.  (95)  and  Eq.  (96) 


are  important.  On  the  other  hand  if  the  current  neutralization 


q z 2 

is  nearly  complete  (e.g.,  £90%  complete)  so  that  uce/“ce  ''<0(e  ), 


then  the  diamagnetic  current  in  Eq.  (95)  also  becomes  important, 


0 2 

whereas  the  terms  proportional  to  (i»>ce)  become  even  less  important. 


To  illustrate  more  clearly  some  of  the  different  types  of  equilibria 
that  can  be  treated  consider  the  following  three  cases: 


0 Z 

Case  i)  - Here  the  current  neutralization  is  poor,  to  /to  MHO 

ce  ce 


% 


and  — v 0(e).  Keeping  only  the  leading  terms  Eqs.  (95)  and  (96)  reduce  to 
n 

e . z 2 e 

ir*  ■ -r  t + 0(6  > <98) 

c to 

ce 


13  , 6 \ Pb  r.b  , n,  2. 

— — (roj  ) t-  V + 0(e  ) 

r 3r  ce  2 z 

c 


(99) 


— b 2 

where  from  Eq.  (55)  one  has  P„  - m yV  +0(e  ).  Since  there  is  no 


contribution  from  the  diamagnetic  current  the  beam  will  be  paramagnetic. 


Case  ii)  - In  this  case  the  current  neutralization  is  approximately 


0 z ^b  2 

50%  complete,  to  /to  v 0(e)  and  — v 0(e  ).  Again  only  keeping  the 

e 


leading  terms  the  field  equations  are 


duo 


e 


ce 


3r 


ce 

z 


ce 


• 2 2 
to  . . to 

vb  * -K  V* 
2 z 2 z 
c c 


+ 0(e3) 


(100) 


13,0. 
7 77  (ra)ce) 


- 2 2 
to  . . (0 

_Eb  yb  + ve 
2 z 2 z 

c c 


+ 0(e2) 


(101) 


Here  again  the  beam  is  paramagnetic.  One  particularly  simple  choice 


~G 

for  V is  just 
z J 


V^(r) 


<i(D  h 

- VbF 

2 z N 


(102) 


pe 


where  ■ constant  and  is  the  fraction  of  current  neutralization. 
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In  general,  however,  F^  may  be  function  of  r. 

Case  iii)  - In  this  case  the  current  neutralization  is  nearly 


complete,  w /w  '0(e~)  and  n,/n  v0(r  ).  Here  one  finds  that 
ce  ee  be 


3w 


ce 


3r 


ce 

z 

ce 


• 2 2 
U)  w 

Vb  + -E£  V 


1 1,1  . p*".  3 An 

e I _ pb  *th 

z|  2 2-  z 3r 

c 2m  yw 


— + 0(t4) 


1 3 , 0 

— t—  (ru>  ) 
r 3r  ce 


e ce 

2 2 
<0  , , (0 

_Pb  vb  + Ve 
2 z 2 Vz 

c c 


+ 0(e3) 


(103) 


(104) 


From  Eq.  (103)  it  is  clear  that  the  beam  may  be  diamagnetic  if  the 

current  neutralization  is  good.  Of  course  if  the  current  neutralization 

0 

is  very  good,  w 'vO  and  the  beam  is  obviously  diamagnetic. 

These  three  cases  illustrate  the  flexibility  of  this  model  in 
treating  a variety  of  beam-plasma  equilibria.  In  the  next  Section 
some  specific  examples  will  be  presented  for  these  three  cases. 

One  other  interesting  case  arises  when  the  beam  is  neither 
paramagnetic  nor  diamagnetic,  implying  that  the  azimuthal  current  is 


zero.  Consequently  is  uniform  and  Eqs.  (103)  and  (104)  reduce  to 


„ 8 , ,r.b„  . , Pth  3lnnb(r) 

0 " “ce(r,V2a-V  - 71=  K 


13  , 8 , wpb 

7-s7  <"V(t)) 2 


2m  y 


v"(1-f„) 

Z N 


(105) 


(106) 


where  it  is  assumed  for  simplicity  that  V^(r)  is  given  by  Eq.  (102) 
with  F^=constant.  It  is  easy  to  show  that  the  other  terms  in  Eqs.  (95) 

and  (96)  can  be  ignored  for  all  values  of  FN  when  J <0(e3',  as  long 

r ||  0 

as  ^ 0(1X  Using  Eq.  (105)  to  eliminate  a)  (r)  in  Eq.  (106),  one 

m c ce 

e 

obtains  an  equation  for  n^(r).  Solving  this  equation  for  n^(r)  one  finds 

\<°> 


2 2 2 
(1+rVr  £) 


(107) 
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which  is  the  density  profile  derived  from  the  Bennett  distribution 
where  the  beam  radius  must  be 

2c  p 


25,26 


rb  = 


th 


to  (0)Vbm  y1/2(1-F  ) 
pb  z e N 


(108) 


If  the  beam  radius  is  not  equal  to  r^,  then  dB^/dr^O  even  though 


n^(r)  is  given  by  Eq.  (107).  B^  is  diamagnetic  when  r^<r^  and 


paramagnetic  when  r^>r^. 


If  E^=0  (i.e.,  there  is  no  return  current),  then  Bz=const.  implies 
that  the  azimuthal  fluid  velocity  of  the  beam  vanishes  so  that  the  beam 
is  nonrotating.  Physically  this  means  that  the  diamagnetic  current  of 
the  beam  exactly  cancels  the  azimuthal  current  of  the  beam  associated 
with  the  electrons  streaming  along  the  helical  field  lines  (plus  any 
azimuthal  particle  drifts,  which  are  generally  negligible).  This 
cancellation  persists  across  the  entire  profile  of  the  beam.  Fig.  1 
illustrates  this  point  by  showing  a typical  beam  electron  trajectory. 

In  particular,  the  sense  of  rotation  of  the  gyromotion  is  opposite  to 
that  of  the  rotation  due  to  the  streaming  along  the  helical  field  lines. 

If  n^(r)  is  given  by  Eq.  (107)  but  FN=const.^0,  it  is  still  possible 
to  have  a nonrotating  beam,  however  the  rotation  of  return  current 
produces  a residual  diamagnetism.  In  this  case  the  radius  of  the  non- 
rotating beam  is  found  to  be 
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Note  that  when  F =0,  r =r  . Finally  if  F„  is  chosen  as  a function  of 
N b b N 

r,  equilibria  can  be  found  with  Bz=constant,  however,  n^(r)  will  not 
correspond  to  the  Bennett  profile  and  in  general  the  beam  will  rotate. 
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From  this  discussion,  it  is  clear  that  a nonrotating  beam  is  a 
rather  special  case  and  that  beams  generally  do  rotate.  Of  course, 
the  physical  consequences  of  this  rotation  in  a stability  analysis 
will  depend  on  the  magnitude  of  the  rotation  frequency  compared  to 
the  time  scale  of  interest. 

IX.  EXAMPLES  OF  TYPICAL  BEAM-PLASMA  EQUILIBRIA 


Two  density  profiles  that  will  be  considered  are  the  Bennett 
profile  and  the  Gaussian  profile.  The  Bennett  profile,  however,  is 
more  diffuse  than  the  Gaussian  profile  (Fig.  2).  This  implies  that 
the  steeper  gradients  at  the  edge  of  the  Gaussian  profile  will  produce 
larger  diamagnetic  currents.  This  will  have  important  consequences 
for  the  beam  rotation  frequency. 

If  Eq.  (102)  is  used  for  V^(r)  where  for  simplicity  FN=const. , 
then  Eqs.  (99),  (101)  and  (104)  are  mathematically  equivalent  and  are 
easily  integrated: 
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where  F =0  for  case  i,  F.^,5  for  case  ii  and  F >.9  for  case  iii  of 
N N 

Sec.  VIII.  For  the  Bennett  profile,  Eq.  (107),  this  can  be  integrated 
to  give 
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For  the  Gaussian  profile  n^r)  is  expressed  as 


nfa(r)  = nb(0)exp{-r2/r^} 


(112) 


and  the  solution  of  Eq.  (110)  for  this  case  is  just 
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0 / -> 
u (r) 
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(113) 


Note  that  in  both  cases  w (r)  increases  linearly  with  r near  the 

ce 

center  of  the  beam  (r<<r,),  and  falls  off  like  r ^ outside  the  beam 

b 

q 

(r>>r_).  co  is  plotted  as  a function  of  r in  Fig.  3,  where  for 
b ce 

comparison  the  density  on  axis,  11^(0),  and  the  beam  radius,  r^,  are  set 
equal  for  the  two  different  profiles.  Although  the  total  current 
is  the  same  for  both  profiles,  the  peak  value  of  Bg  is  ^25%  higher 
for  the  Gaussian  profile  than  for  the  Bennett  profile  simply  because 
the  Bennett  profile  is  more  diffuse.  By  choosing  appropriately, 
the  curves  will  correspond  to  any  one  of  three  cases  discussed 
in  Sec.  VIII.  With  certain  values  for  the  parameters  this  model  can 
properly  describe  self  fisids  as  large  as  ^10%  of  the  applied  field 
for  a Gaussian  profile  (Fig.  4).  B.  could  become  even  larger  for  a 
density  profile  with  a sharper  boundary  than  the  Gaussian  profile. 
Turning  to  the  case  where  J0=O,  it  is  interesting  to  examine 

g 

how  the  axial  self  magnetic  field,  B^ , and  the  beam  rotation  frequency, 
wR,  are  modified  by  varying  r^.  Here  the  beam  rotation  frequency,  wR, 
is  defined  as 
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If  r =r,  then  B will  vanish  for  all  r,  and  if  r.=r.  then  will 
b b z ’ b b R 

g 

vanish  for  all  r.  B (r)  and  o>  (r)  are  plotted  in  Figs.  (5)  and  (6) 
for  four  different  values  of  r^  where  the  beam  density  follows  a 
Bennett  profile.  Observe  that  the  diamagnetic  beams  are  thinner  and 
generally  rotate  faster  near  the  axis  than  the  paramagnetic  beams. 
Although  Fn=0.5,  the  diamagnetic  beams  considered  here  differ  from 
the  paramagnetic  beams  of  case  ii  in  Sec.  VIII  because  the  beam 
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radius  is  smaller.  Figs.  (7)  and  (8)  display  the  corresponding  curves 

g 

for  Bz(r)  and  wR(r)  when  n^(r)  follows  a Gaussian  profile.  Note 
that  (^(r)  depends  critically  on  the  shape  of  the  density  profile  while 
only  has  a weak  dependence  on  n^(r).  In  particular,  the  beam 
rotates  in  the  same  direction  across  the  entire  profile  of  the  beam 
for  the  Bennett  case  for  any  r^,  whereas  the  beam  can  change  its 
direction  of  rotation  at  the  outer  edge  of  the  beam  for  certain  values 
of  r^  when  the  density  profile  is  Gaussian.  Although  the  behavior  of 


B (r)  about  r =r  is  interesting  its  magnitude  is  small  compared  to 

Z D D ^ 

the  applied  field  since  it  was  assumed  that  w9  /wZ  <0(e)  and  — — < 0(e) 
r ce  ce'v  m c 'v 

6 z 6 

Larger  values  of  p^,  /m  c or  u>  /w  will  of  course  make  these  effects 
th  e ce  ce 


more  important. 

If  the  current  neutralization  is  better  in  the  interior  of  the 
beam  than  at  the  edge  of  the  beam,  then  a larger  portion  of  the  net 
current  flows  on  the  surface  of  the  beam.  This  situation  can  arise 


when  r,  co  /c>l  and  the  return  current  has  not  yet  diffused  out  of 
b pe  J 
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the  beam  channel.  * This  situation  is  easily  modeled  by  allowing 


Fn  to  be  a function  of  r.  One  convenient  choice  is 


Vr)  - Vb<r>/nb<0) 


Figs.  (9)  and  (10)  show  the  effect  that  this  has  on  B (r)  for  the 

u 


Bennett  profile  and  the  Gaussian  profile  respectively  for  different 


values  of  F^.  When  F^=l,  B0  still  reaches  an  appreciable  peak  value 


since  the  current  is  only  completely  neutralized  at  the  center  of  the 
beam.  Furthermore  with  the  current  flowing  nearer  to  the  edge  of  the 


beam,  B.  peaks  further  from  the  center  of  the  beam  than  if  F.,(r) 
U N 


( 


were  constant. 


Although  the  choice  of  FN(r)  given  in  Eq.  (115)  is  instructive, 
any  reasonable  function  of  r can  be  used  for  F^(r)  in  order  to  fit 
any  appropriate  profile  for  Bg(r).  Thus  both  n^(r)  an^  BQ(r)  can  be 
chosen  independently.  This  fact  makes  this  Vlasov  treatment  a very 
powerful  method  for  modeling  beam-plasma  systems. 

X.  CONCLUSIONS 


The  Vlasov  equilibrium  model  that  has  been  developed  in  this 

paper  can  describe  a very  broad  class  of  warm  beam-plasma  systems  in 

an  applied  field,  B^;  but  the  entire  treatment  depends  crucially  on 

the  smallness  of  the  electron  gyroradius  compared  to  the  beam  radius. 

This  condition,  however,  is  easily  satisfied  in  typical  systems  where 

B^>1  kg  and  r.  >1  cm.  Assuming  that  this  condition  on  the  gyroradius 

is  satisfied,  the  procedure  outlined  in  this  paper  can  be  used  to  model 

any  charge  neutralized  beam-plasma  system  with  reasonable  density 

and  azimuthal  self  magnetic  field  profiles.  The  density  profile  and 

azimuthal  field  profile  can  be  chosen  independently  since  the  theta 

component  of  the  V*Jjl  equation  can  be  satisfied  by  adjusting  the  fractional 

current  neutralizational  factor,  F„(r) . In  other  words  the  (VxB)Q 

equation  can  be  used  to  determine  FN(r)  given  the  functional 

dependence  of  n^  and  B0  on  r.  Of  course,  if  F^fr)  is  known,  then 

the  equation  can  be  used  to  find  B (r). 

By  making  the  simple  choice  of  FN(r)=const.  and  using  the 

equation  to  solve  for  Bn(r)  and  B (r) , it  was  possible  to  investigate 

o z 

the  dependence  of  the  beam  rotation  frequency,  ojR(r) , and  the  axial 

g 

self  magnetic  field,  Bz(r),  on  the  density  profile  of  the  beam,  n^(r), 

and  the  fraction  of  current  neutralization,  F„.  When  F >0.5  it  was 

N N 
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g 

found  that  B^  can  be  diamagnetic  for  a sufficiently  warm  beam  but  when 

g 

F <0.5,  B was  paramagnetic.  In  general  diamagnetic  beams  were  also 
N Z 

found  to  be  thinner  than  paramagnetic  beams  because  of  the  increased 
pinching  force  when  J0<0.  It  was  discovered  that  depends  strongly 

on  the  shape  of  the  density  profile  and  that  the  Bennett  profile  with 

<\, 

radius  r,  (Eq.  109)  corresponds  to  the  case  of  a nonrotating  beam,  i.e., 
b 

iuD(r)=0  for  all  r.  However,  this  nonrotating  beam  is  only  one  special 
K 

case  of  the  general  type  of  beam-plasma  systems  described  by  this 
equilibrium  model. 

Because  of  the  flexibility  in  the  choice  of  n^r)  and  B((r) 
mentioned  above,  this  equilibrium  model  will  be  very  valuable  for 
studying  the  effects  of  self  magnetic  fields  and  radial  inhomogeneities 
on  the  stability  of  beam  plasma  system.  Since  the  model  also  allows 
for  a thermal  spread  in  momentum,  thermal  effects  on  the  stability  of 
such  systems  can  also  be  studied.  Furthermore,  as  mentioned  in 
Sec.  Ill,  the  beam  electron  orbits  found  in  Eqs.  (41),  (42)  and  (43) 
and  background  electron  orbits  found  in  Eqs.  (47),  (48)  and  (49) 
can  be  utilized  in  calculating  the  perturbed  distribution  functions 
for  a stability  analysis  using  the  method  of  characteristics.  This 
analysis  will  be  the  subject  of  a future  paper. 


Fig.  1 Trajectory  of  a typical  beam  electron.  The  beam  electron 

gyrates  about  ^ in  the  opposite  sense  that  it  rotates  about 
the  z-axis  while  streaming  along  B. 


Fig.  2 Beam  density,  n^(r),  as  a function  of  the  radial  position,  r. 

Both  Bennett  and  Gaussian  profiles  are  shown,  with  u^^(0) 

and  r,  identical  for  both  cases, 
b 


Fig.  3 Azimuthal  self-magnetic  field,  Bg(r),  as  a function  of  the 

radial  position,  r.  Here  the  fractional  current  neutralization, 

Fn,  is  a constant  (0<FN<1).  Examples  of  the  three  cases 

discussed  in  Sec.  VIII  are  obtained  when  F <0.1  (case  i) , 

N'v 

F,/\0.5  (case  ii)  and  F >0.9  (case  iii) . The  parameters  of 

N N''' 

the  system  are  identical  for  both  the  Bennett  case  and  the 

Gaussian  case:  r =1  cm,  w =9.9xl0^rad/sec;  wz  (°°)  = 3.0xl0^rad/sec; 
b pe  ce 

n^(0) /np(0)=0.01;  y=2;  pt^/mec=0.1. 


Fig.  4 Example  of  a case  where  the  azimuthal  self  magnetic  field 
is  greater  than  10%  of  the  applied  field.  Here  a Gaussian 
density  profile  is  used  and  r^=2.0  cm;  w e=6. 75x10^rad/sec; 

7 10 

a>ce (°°)=3.0xl0  rad/sec;  F^=0;  n^(0) /np(0)=0.04;  y=2.0; 

Pti/V"0-1- 

Fig.  5 Axial  self-magnetic  field,  B^(r),  as  a function  of  radial 

position,  r,  for  a Bennett  density  profile  with  various 

values  of  the  beam  radius,  r,  . Note  Bs(r)=0  for  r =r.=1.0  cm. 

b z b d 

The  other  parameters  of  the  system  are  the  same  as  in  Fig.  3 


with  Fn=0.5. 
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Fig.  6 


Beam  rotation  frequency,  w (r),  as  a function  of  radial  position, 

K 

r,  for  a Bennett  density  profile  with  various  values  of 
the  beam  radius,  r,  . Note  u)„(r)=0  for  r =r  =0.7  cm. 

D K b D 

The  other  parameters  of  the  system  are  the  same  as  in  Fig.  3 
with  F =0.5. 
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Fig.  7 


Fig.  8 


Axial  self-magnetic  field,  B^(r),  as  a function  of  radial 
position,  R,  for  the  case  of  a Gaussian  density  profile  for 

g 

various  values  of  the  beam  radius  r,  . Here  B (r)  is  nonzero 

□ z 

for  all  r^.  The  other  parameters  of  the  system  are  the  same 
as  in  Fig.  3 with  FN=0.5. 


Beam  rotation  frequency,  u_(r),  as  a function  of  radial 
position,  r,  for  a Gaussian  density  profile  with  various  values 
of  the  beam  radius,  R^.  Note:  there  is  no  nonrotating  case 
for  any  R^.  The  other  parameters  of  the  system  are  the  same 
as  in  Fig.  3 with  FN=0.5. 


Fig.  9 Azimuthal  self  magnetic  field,  B (r) , as  a function  of  radial 

V7 

position,  R,  for  a Bennett  density  profile.  Here  F„(r)= 

N 

FNn^(r)/n^(0)  for  ?N  =0.0,  0.5  and  1.0.  The  other  parameters 
of  the  system  are  the  same  as  in  Fig.  4. 


Fig.  10  Azimuthal  self  magnetic  field,  B (r) , as  a function  of 

radial  position,  R,  for  a Gaussian  density  profile.  Here  FN(r)< 
F^n^(r) /n^(0)  for  FN=0.0,  0.5  and  1.0.  The  other  parameters 
of  the  system  are  the  same  as  in  Fig.  4. 
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